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Abstract 

A semiclassical model is used to investigate the possibility of selectively exciting one of two 
closely spaced, uncoupled Raman transitions. The duration of the intense pump pulse that creates 
the Raman coherence is shorter than the vibrational period of a molecule (impulsive regime of 
interaction). Pulse shapes are found that provide either enhancement or suppression of particular 
vibrational excitations. 

PACS numbers: 42.65.Dr 42.50.Hz 33.20.Fb 82.53.Kp 
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I. INTRODUCTION 



Studies of molecular dynamics accompanying excitation, ionization, or dissociation give 
rise to questions concerning the control of specific motion in molecules. Control using ul- 
trafast lasers and adaptive learning algorithms, e.g., genetic- type algorithms , allows 
one to optimize feedback and achieve predetermined goals. Series of successful experiments 
implementing this technique offer proof that it is possible to exercise quantum coherent con- 
trol over a range of processes. Two such experiments employ stimulated Raman scattering 
in liquid methanol [2] and gas phase carbon dioxide JjJ . In the former case a liquid medium 
is excited with a 100 fs pulse, which is equivalent to non-impulsive Raman excitation and, 
therefore, makes propagation effects important. In the latter case the selective Raman exci- 
tation of vibrational modes is achieved with a pulse of duration within a typical molecular 
vibrational period. In this case the seed for the Stokes transitions is within the pump pulse. 
In the experiment with methanol the problem of intra-molecular couplings of vibrational 
modes is discussed using a comparative analysis with experiments on benzene and deuter- 
ated benzene molecules. The controlling optical field may be constructed using coherent 
quantum control theory, see, e.g. 0, 0|. An alternative approach involves a search for an 
analytical pulse function. In [aQ] an antisymmetric phase function is proposed, which in- 
verts the sign of the electric field at a given energy, inducing constructive interference of the 
off-resonant components of the spectrum and therefore maximizing the transition amplitude. 
Selectivity of the two Raman levels of pyridine, both of which lie within the bandwidth of 
the excitation pulses is achieved in CARS signals by positioning the phase step at each of 
the two peaks jj]]. In [l| the use of adaptive techniques in FAST CARS is explored for an 
identification of related biological objects such as bacterial spores. 

The Raman excitation of a single vibrational mode having frequency u normally requires 
a pulse with lot < 1, where r is the pulse duration. Optimal excitation may occur if the 
intensity is modulated at a frequency that corresponds to the transition frequency. For any 
two vibrational modes uix, U2 within the pulse bandwidth, {u\ — u)2)r <C 1, the spectral 
width is too large to resolve the lines. In the present paper we propose a method for 
achieving selective excitation using broadband pulses. In the frequency domain we introduce 
an intensity envelope that vanishes for a modulation frequency of the pulse equal to the 
frequency of vibration that we would like to suppress. The Fourier transformed field, when 
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applied to a molecular system, provides negligible excitation of that particular frequency 
together with significant excitation of another frequency. This picture is valid for weak 
fields. In strong fields the effect is not obvious and numerical analysis reveals that, for some 
intensities and the designed pulse shape, it is still possible to optimize one transition versus 
another. 

The paper is organized as follows. In the second section the problem is formulated based 
on a model of two, two level systems. An analytical function for the intensity envelope is 
proposed and the equations of motion for the state amplitudes are obtained. The third 
section contains results of numerical calculations that determine the field parameters that 
lead to selective excitation. The paper is ended with a short summary. 

II. BASIC FORMALISM 

In our model a molecule is described by two, two-level quantum systems, each representing 
a normal Raman-active vibrational mode. The molecular medium is represented by an 
ensemble of such systems with no relaxation processes taken into account. Each two-level 
system interacts with an intense off-resonant femtosecond pulse that initiates stimulated 
Raman scattering via an off-resonant interaction with a virtual state. The duration of this 
pump pulse is shorter than a typical vibrational period of a molecule. In this case the 
frequencies of both two-level systems are within the bandwidth of the pulse and the Stokes 
component of the field is supplied by the same pulse. The coherently excited ensemble is 
analyzed by a weak probe pulse, not considered in this work, applied after a time delay 
shorter than the coherence time. The goal of the present work is to determine a pulse shape 
that provides selectivity for the excitation of one of the two-level systems. 

A semiclassical model of laser-molecule interactions is used. The model is represented 
schematically in Fig.l where u>2i is the frequency of the first two- level system and CJ43 that 
of the second system. Initially only lower levels |1 > and |3 > of both two- level systems 
are populated, and the populations of these levels are equal. The time evolution of two, 
two-level systems is described in terms of probability amplitudes, which are written in the 
interaction representation 
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Here cxjiv is the frequency of a single level, such that, e.g., (a 2 — oli)uj = U21, is a Rabi 
frequency, fijb is a dipole moment matrix element, E p0 (t) is the pulse envelop, A is the 
detuning of the frequency of the pulse from the frequency of the virtual state \b >. Note, 
that the pulse envelope E p0 (t) is the same for all transitions. The Rabi frequencies may 
differ owing to different dipole moment matrix elements. 

The system of coupled differential equations (1) is derived from the time-dependent 
Schrodinger equation with Hamiltonian: 
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where oj p is the laser field carrier frequency. By adiabatically eliminating state \b > within 
the rotating wave approximation, we arrive at Eqs (1). In this work we discuss the case of un- 
coupled two-level systems such that the probability for the population flow from one system 
to another via the external field is zero. Then Eqs.(l) are represented by two independent 
systems of coupled differential equations with two variables. 

Coherent excitation of a molecular medium induces a vibrational coherence | P12 1 or | /O34 1 . 
When a probe field is applied on the 1-b transition, the coherence |pi 2 | serves as a source 
for the generation of a field on the 2-b transition. Thus it is of prime interest to calculate 
I />i2 1 and I /O34 1 . The goal of this paper is to choose a pump pulse such that | P12 1 is suppressed 
while I P34 1 is enhanced. 

We propose an analytical function for the intensity envelope which is included in the 
dynamical equations (1) for the probability amplitudes. It is easiest to choose this function 
in the frequency domain. To suppress excitation at frequency cu 2 i and enhance excitation at 
frequency u> 43 , we choose 

/» = I e~^ «) 2T2 (l - e ^(— ) 2 ^ 2 ) , (3) 

where T and Ti are free parameters. When the modulation frequency of the pulse u is equal 
to W43 the intensity approaches its maximum I for a sufficiently large parameter Ti, 
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I(oj 43 ) = J (l - e~ A ^ 2 ), Alo = lo 43 - U21. (4) 

For iv equal to 1V21 the intensity is zero, I(iV2i) — 0. The intensity envelope as a function of 
frequency is drawn in Fig. 2. The frequencies of the two-level systems are iv 2 i—l and cj 43 =1.1 
in frequency units of a> 2 i- The intensity of the field at iv i3 and frequency region near iv 2 i over 
which I((v) is approximately zero depend on T x . The larger Ti, the greater is the selectivity 
for suppressing the oj 2 i frequency. 

The inverse Fourier transform of the spectral density (3) is a complex function. To arrive 
at a physically acceptable temporal pulse function, we take the real part of the inverse 
Fourier transform, given by 

lit) = J ((V2T)- 1 e-&cos(u 43 t) - (v^r)" V*" 2 ^ 1 -^^ cos((cv 21 - Alo^)^ , (5) 

where r 2 = T 2 + T 1 2 . The expression for the field (5) is inserted in Eqs.(l) for the calculation 
of the probability amplitudes. 

The solution of Eqs.(l) in the limit of a weak field using perturbation theory is 

«4 = ^ I(t)e^dt = i^lM = <^J (1 - e"^ 2 ), 

(6) 

a 2 = Hfie^dt ~ e ~ T2 ^ + ^ 2 + e ~ A ^ T2 e ~ r2 ^ 2 w 0. 

The Fourier transform represented in Eq.(5) is not identical to Eq.(3), since we took the 
real part of the Fourier transform to arrive at (5). It now contains "counter-rotating" 
terms, which are small for the chosen pulse shape. Thus, by construction, we have totally 
suppressed the 1-2 transition in the weak field limit. On the other hand, the excitation of 
the 3-4 transition is still weak owing to the perturbative nature of the solution. 

Polyatomic molecules often possess several or many Raman active modes with frequencies 
close enough to be within the bandwidth of the pulse. In order to enhance a single vibration 
and suppress other vibrations the function for the pulse may be constructed as a product of 
several terms 

/» = J e-( w - w ^ 2T % (l - e -(<"-«; W) . (7) 
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III. NUMERICAL RESULTS 



In this section we discuss the results of numerical calculations based on the exact solution 
of Eqs.(l). The numerical studies reveal the influence of the field parameters on the efficiency 
of the excitation of the two-level systems. Parameters for the system are taken from the 
experimental data on impulsive excitation of vibrational modes in the molecular gas CO2 
P|. In CO2 the frequencies of two selectively excited Raman modes are 36.8 and 42 THz. 
The FWHM of the applied intense pulse is taken equal to 18 THz. In our calculations the 
frequency uj 2 i is set equal to unity; in these units the frequency U43 is equal to 1.1. From 
experimental data, we estimate that the parameter T is about equal to 3, in frequency units 
of u>2i- The intensity of the field is determined by the parameter I . The parameter Ti is 
related to the width of the spectral dip in I{ui) centered at frequency U21 ■ Although a value 
of 3> T would provide the best selectivity, the choice for the parameter T\ is strongly 
restricted by the requirement that the duration of the applied pulse be within a typical 
molecular vibrational period. It turns out that even for such values of 7\, it is possible to 
selectively excite one transition. 

We calculated the population distribution and the coherence of the excited and suppressed 
two-level systems as a function of T\ as shown in Fig. 3. Bold solid and dotted lines depict 
the absolute value of the coherence |p 34 | of the excited system and |pi 2 | of the suppressed 
system. Populations of the upper levels of both systems are shown by bold dashed and 
dot-dashed lines. The intensity of the field Iq is 7r/8 which corresponds to a weak field in 
our calculations (but not to the perturbative regime). For the value 7\ = 13 the population 
of levels of the 3-4 system is 0.25, and the coherence is optimal for the given intensity of the 
field and dipole moments, (fii = 1). The duration of the laser pulse corresponding to this 
value of the parameter 7\ is about 200 fs which does not satisfy the necessary requirements 
on pulse duration. According to Fig. 3 for smaller values of the parameter T 1; corresponding 
to shorter pulses, the coherence of the 3-4 system is significantly reduced with a simultaneous 
increase of the coherence of the 1-2 system. Optimal values for the coherence of both the 
1-2 and 3-4 systems were found for T\ <T through a search over different intensities of the 
field. 

In Fig. 4 the coherence is plotted as a function of the intensity of the field for parameters 
T=3 and T\ = 3. Coherence | P34 1 of the 3-4 system is represented by a bold solid line 
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and coherence | yf>i2 1 of the 1-2 system by a bold dashed line; thin lines show populations of 
the upper levels of both two-level systems. For the intense fields, coherence of the excited 
and suppressed systems possess somewhat chaotic structure. Several values of the intensity, 
e.g. Iq = 7r and I = 1.757T give rather low coherence of the 3-4 system but maximum 
coherence of the 1-2 system. A desired solution for maximum coherence of the excited 3-4 
system is achieved for the intensity coefficient Iq = 2.087T. This is the result of redistribution 
of population within that two-level system: half of population is transferred to the upper 
level providing maximum coherence. The corresponding coherence of the 1-2 system at this 
intensity is nearly zero, where most population remains in the lower level. For the intensity 
I = lAir the picture is similar, however the low coherence of the 1-2 system is due to nearly 
complete population transfer to the upper level. 

The goal of control of the coherence of two uncoupled two-level systems is achieved with 
a pulse shape possessing a broad spectral dip at the suppressed frequency and a suitably 
chosen intensity of the field. This technique allows one to use pulses of duration T to 
selectively excite transitions having frequency separations Au < 1/T. Our results should 
not be taken to imply that one can spectroscopically determine frequencies to better than 
the inverse temporal width of the pulse. On the other hand, if the frequencies are known 
from previous measurements, it is possible to suppress one transition and enhance the other 
by the method outlined above. 

Had we taken a frequency profile centered at 0^43 with T = 3, the curves for the coherence 
and populations of the excited and suppressed systems would differ qualitatively from those 
show in Fig. 4. The desired selectivity could not be achieved. 

The time-dependence of the coherence, populations and the field is shown in Fig. 5 for T = 
3, Ti = 3, Iq = 1.47T. The pulse duration is about 50 fs (thin solid line). It induces oscillations 
in the population distribution (thin lines) which lead to oscillations of the coherence of the 
two- level systems (bold lines). At long times the coherence and populations of levels achieve 
stationary values. 

We have carried out some preliminary calculations, including the coupling between vibra- 
tional modes via an external field used as a control mechanism. If the direct excitation of a 
particular vibrational mode is weak owing to a weak oscillator strength, the excitation may 
be enhanced through the coupling to the other Raman active vibrational modes. Details 
will be published in 
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IV. SUMMARY 



We have presented a semi-classical description of stimulated Raman scattering involving 
the selective excitation of one of two closely spaced vibrational modes. The dynamics is 
described in terms of probability amplitudes and depends on the shape of a femto-second 
laser pulse where duration is shorter than a typical molecular vibrational period. We propose 
an analytical function for the shape of the intensity envelop of the pulse that allows for the 
selective excitation of a predetermined vibrational motion with simultaneous suppression 
of an unfavorable one. The pulse leads to maximum coherence of a desired vibrational 
transition and consequently to maximum gain into Raman side bands when a probe pulse 
is applied. This pulse function may be used as an initial guess for the control of bond 
excitation in chemistry as well as a fit function within an adaptive learning algorithm. 
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Ib> 




FIG. 1: Schematic picture of a model system consisting of two, two-level systems having frequen- 
cies LJ21 and Initially, the lower levels are populated evenly. The uncoupled transitions are 
driven by an off-resonant femtosecond pulse. 
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FIG. 2: Intensity spectral profile as a function of frequency for T\ = 10, 5, 3 (dashed, dotted, and 
solid lines). In the insert the intensity envelop as a function of time is presented for T\ = 10,3. 
All frequencies are in units of a>2i and times in units of uj^i- 
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FIG. 4: Intensity dependence of the coherences \p±2\ and \pu\ and upper states populations P22 
and P44. Maximum coherence of the 3-4 system and negligibly small coherence of the 1-2 system 
are observed for 1$ = 1.4 and 2.08 in the intensity region shown with T = 3, T\ = 3. 
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